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to a considerable delay of its publication, I have decided now to publish it in a preprint form. 
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1. Introduction 

Since a long time it has been well known that the number of Newtonian differential 
equations of motion can be diminished by making use of the existence of some me¬ 
chanical conservation laws. In the middle of the 19 th century a new problem of this 
kind was posed. If the original set of equations of motion are the Euler-Lagrange 
equations of a Lagrangian, and the dynamical system admits conservation laws that 
can be used to reduce the number of these equations, is then always possible to find a 
new Lagrangian such that the reduced system of differential equations can be derived 
as Euler-Lagrange equations of the new Lagrangian? 

As is known, cf. T], the first solution to the problem was given in 1876 by 
E. J. Routh who showed that when the corresponding conservation laws resulted 
from the occurrence of cyclic dynamical variables in the original Lagrangian, it was 
the Routh function that was the Lagrangian for the reduced system. If, however, 
the law responsible for the reduction of the system is the conservation of energy, the 
Routh method applied directly to the original Lagrangian does not work at all. This 
case required a separate treatment that was given in 1886 in a book by K. G. J. Ja¬ 
cobi, where a variational principle leading to differential equations satisfied by spatial 
trajectories in the configuration space of the dynamical system was formulated under 
the assumption that the original Lagrangian did not depend explicitly on the time 
variable. The proof of Jacobi was based on the Maupertuis principle of least action. 
This fact may be one of the sources of a terminological confusion which appears in 
many contemporary text-books on analytical dynamics, where the Jacobi principle is 
named Maupertuis principle, despite the fact that the latter is, for holonomic dynam¬ 
ical systems and for £ / 0, equivalent to the Lagrange equations of the second kind 
which determine the motion of the system, whereas the Jacobi principle determines 
only the orbits of the motion. 

In 1994, in J2], the present author together with P. Jaranowski have posed and 
solved, as they have named it, the inverse Jacobi problem: under which conditions 
imposed, can one restore the original motion when starting from a variational principle 
leading to orbits? 

In |2] an attempt was also made to derive ab initio the standard Jacobi principle 
without making any use of the Maupertuis principle. During my seminar talks on the 
results obtained in |2J, I realized that the “new” derivation of the Jacobi principle 
was rather complicated to convey it to the audience. In 2001, I found a very simple 
derivation of the Jacobi principle, published later in 2 ]. It makes use of the Routh 
method applied to an in a suitable way transformed Hamilton’s action. It is so simple 
that it must have been undoubtedly known to some people before, although I could 
not find any references to it. From the point of view of methodology, this derivation 
is more suitable for a classroom than the traditional one, because it solves two akin 
problems in the same way. 

The main objective of the article is to demonstrate that there exist mappings be¬ 
tween dynamics of classical systems on one side and variational principles for geodesic 
lines in stationary spacetimes of general relativity on the other. The construction of 
the mappings is based on classical Routh’s and Jacobi’s reduction procedures and 
on corresponding inverse procedures which are proposed by the present author and 
P. Jaranowski in [2] and 0 . All these procedures are general theoretical methods 
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that belong to analytical dynamics. A review of them is presented in sections 2, 3, 
4, and 5, mainly in order to fix the framework which will be employed in the next 
sections. 

Sections 6, 7, and 8 present the classical Jacobi procedure in the working. In 
Sec. 6, the relation between two widely known actions for geodesics on manifolds is 
interpreted in terms of the Jacobi reduction of the “quadratic” action into the other 
one. Section 7 repeats the elementary text-book example of the Jacobi reduction of a 
Newtonian, holonomic Lagrangian into the Lagrangian describing orbits as geodesics 
in the kinetic energy metric. The action for geodesic lines in a stationary Lorentzian 
manifold in the coordinate time parametrization is in Sec. 8 Jacobi reduced into an 
action defined on the constant time hypersurface. All the examples considered in 
these three sections are at the end of Sec. 8 reinterpreted in terms of mappings of 
some of the dynamics into the other ones, and the equivalence of some of the dynamics 
is exhibited there. 

The inverse Jacobi procedure is applied to the action considered in the previous 
section in Sec. 9. Its result is an action determining affinely parametrized geodesics, 
and the metric coefficients in this action are time independent. This fact enables 
one to form a composition of the Routh reduction with the inverse Jacobi procedure 
performed just at the beginning of the section. The outcome is a dynamics which 
is equivalent to the Newtonian dynamics considered in Sec. 7. This enables one to 
continue the discussion led at the end of the previous section and to construct two 
closed loops of mappings that alternatively can be considered as a commuting diagram 
of mappings between all the dynamics dealt with in this article. Two of the branches 
in this diagram may be regarded as generalizations of the correspondences between 
dynamics that were already discussed in the literature, cf. and 0, but by methods 
that are rather particular, and without any reference to general principles of analytical 
dynamics. 

In the text which follows, an abbreviated notation is used, in accordance with 
which expressions like e.g. (q l ,t stand for sequences (g 1 , q 2 ,..., q n , q 1 , ..., q n ) or, 
depending on ranges in which the indices vary, for some other sequences of a similar 
type. The summation convention is employed throughout the article. 

2. Routh’s theorem 

Let 

t2 

W[g a ]=/ C(q l (t),q l3 {t),t) dt (1) 

ti 

be an action functional describing a dynamical system in a configuration space Q ra+1 . 
The local coordinates g“, where a,/3 = 0, 1, of a point in Q ra+1 are functions of 
time, q a = g“(i), called the motion of the system in (Q” +1 . Let further the Lagrangian 
C be non-degenerate. The form of the action (1) was written down in accordance with 
the assumption that dC/dq° = 0, i.e. with the fact that the variable q° is a cyclic 
one. 

From this assumption it follows that 
BC, 

Po = := Vq (g , (*),9°(*),? J (i),i) = const, (2) 


where i.j = 1, 2,..., n. 
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Then, cf. JI], 

1. Equation 0 can be solved with respect to the variable q° leaving us with a 
relation of the form 

q°(t) =(j)[po,q l (t),q J (t) 1 t) , 

where po is an arbitrary, but fixed, value of the integration constant. As a result, 
the variables (q°(t), q°(7)) can be eliminated from the system of the n +1 original 
Lagrange equations. 

2. The remaining n differential equations for the variables q l (po, t) are again Euler- 
Lagrange equations of an action integral 

t 2 

Wpo [<f ] = J L po ( q l (t ), q 3 (t ), t) dt, (3) 

tl 

where L po is defined as 

L P0 ( q\q 3 ,t ) = U (q\4’(p 0 ,q k ,q\t),q : ’,t) , (4) 

and where TZ denotes the Routh function 

^(<7 *,<f A\t) := C {q\q°A\t) - q°p 0 . 

3. After the Euler-Lagrange equations corresponding to the action 0 have been 
solved for q l (po, t), one can find the function q°(po , t) by solving the differential 
equation 

<7° =(5) 
where the function </>(po, t) is a solution of the equation 

Vo (V (po, t ), </>(p 0 , t ), q J (p 0 , t ), = po (6) 

into which the now known functions q l (po,t) and q J (po,t) are substituted. 


3. The Routh inverse procedure 

In order to determine the complete motion described by (q°,q l ), the knowledge of a 
pair of functions (L po ,V o), and of a constant po was necessary. A natural question 
now arises whether this information is also sufficient to determine the functional form 
of the original Lagrangian C provided the triple (L po ,Vo,Po) is known. 

The answer to the question just posed is positive, and the proof proceeds as follows. 

1. Suppose that a function Vo is given. Any Lagrangian C(q z ,q°,qi,t) such that 


dC 

dq° 


Vo {q\q°,q J ,t) 


(7) 
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is of the form 

+ A(q l ,q\t), (8) 

where 

IWA 0 A 3 ,t) = J V 0 {q\q°,q J ,t)dq°, 

and A is a quite arbitrary function of the arguments shown in ©■ 

2. The arbitrariness of A is removed by the requirement that the Routh procedure, 
which starts from the assumption 

Vo {q l (t),q°{t),q J (t),t) =p 0 = const, (9) 

if applied to 0 , lead to the now known Lagrangian L po (q l , q\t). As a result, 
one obtains 

A {q\q J ,t) = L po (q\q J ) - I (q l , p(p 0 , q k , q\t) q J ,t) + <p(p 0 , q l , q\t)p 0 , (10) 

where the function ip is defined in an implicit way by the equation 

Vo {q l {t),ip,q>{t),t) =p 0 . (11) 

Of course, the value of the parameter po in Eqs. m and m must agree with 

that entering the known Lagrangian L po (q l , q J ). 

3. The final functional form of £, obtained in consequence of substituting Eq. m 
into Eq. ©, is 

9° 

V(.q\q° 1 q\t) = L po {q\q\t)+^(po,q\q J ,t)po+ j V 0 (q\ n,q J ,t) dn. 

v(po,q\q j J) 

( 12 ) 

Depending on the number of solutions for ip admitted by Eq. m. the solution 
0 of the inverse problem may not be a unique one. 


4. The Jacobi principle 

Let us consider now an action functional of the form 

*2 

W[q]= J L(q\t),q°{t))dt. (13) 

tl 


The form above is equivalent to = 0 which implies the energy conservation law 
G{q\q j ) =E, where 


r)T 

G(M) = ? W -L 


is the energy function, and E is the energy constant. 


( 14 ) 
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In order to bring the action to a form to which the Routh formalism may be 
applied, a transformation of the parameter: t —> r is performed, defined as t = 9(r ), 
where O' (t) ^ 0, and 8 is a meanwhile unknown function. The action m transforms 
then into 

T2 


W[0, x l ] = J A (x l 

Tl 

where 

A = 

and 

x l (r) := 
x'\t) := 


x n (t)) dr , (15) 

i ( X<(T)> ^) 0,{T) ’ (16) 

9 l (%)), ( 17 ) 

g l (0(r))0'(r). (18) 


The new Lagrangian A determines a system of n +1 degrees of freedom described 
by n + 1 independent variables ( 9,x l ) being functions of a parameter r. (Notation 
like a/ = ^ etc is applied here). 

The Lagrangian A is a homogeneous function of degree one in the variables (9 ', a;'*). 
The appropriate variational principle determines thus only n independent differential 
equations of motion regardless of the fact that the system is described by n + 1 
dynamical variables. The Lagrangian A does not explicitly depend on 9. Therefore, 
this variable plays here the same role as q° did in the case of the Lagrangian C 
discussed before. Equation 0 reads now 


Po 


dA 

~d9' 


L 




8'{t) dq k y ’ 



G 



x ,j (t) 

Hr) 


■ ( 19 ) 


Thus Vo = —G(x\ ), and p 0 


—E. Therefore, we have to solve the equation 


G 



= E 


( 20 ) 


with respect to 0', prior to starting with the Routh formalism. 
Writing the solution as 


9'{t) = (j) E (x 1 (t),x ,j (t)^ , (21) 

we are prepared to transform A to a corresponding Routh function which is denoted 
now by Le, 
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L E {x\x' j ) 


A (x l ,(t> E (x J ,x' k ),x' 1 ^ - p 0 (j)E{x\x' J ) 


L x l 


4> E (x k ,x r ) 


■E 


4> E {x r ,x' ) 


( x k X ' 1 \ 

dq l y ’ (t> E {x r ,x' s ) J 


( 22 ) 

(23) 


The Lagrangian L E , for the first time derived by Jacobi, describes a reduced 
dynamical system which resulted from eliminating the information about the time 
evolution from the original system with the Lagrangian L. In other words, the vari¬ 
ables q l (t) which enter L, after the corresponding equations of motion are solved, 
describe motions of the system in Q™ which are curves in Q" parametrized by the 
Newtonian time t. On the other hand, the variables x 1 that enter L E describe tra¬ 
jectories (i.e. spatial paths) of the system; these trajectories are only loci of points 
in Q 71 . As far the computations that determine the form of the Lagrangian L E are 
concerned, the expression (l22l) is, in my opinion, more suitable for practical compu¬ 
tations than the usually quoted expression J2B. It is worthwhile to note that the 
original Lagrangian L provides information about the form of its energy function 
G, whereas this piece of information is lost from the reduced Lagrangian L E ; from 
Eq. Ill 41) it follows that its “energy” function identically vanishes, i.e. no energy - no 
time evolution. 

One can show, cf. f2j, that objects introduced in this section have the following 
properties. 

1. The function cf> E is homogeneous of degree one in the variables x' 1 , which means 
that the relation m is covariant with respect to reparametrizations r —> t' . 

2. This in turn implies that also the Jacobi Lagrangian L E is a homogeneous 
function of degree one in the variables x' 1 . 

3. The rank of the Hesse matrix of L E is equal to n — 1. 

Points 2 and 3 mean that the Lagrange equations 

__ dL E _ d_ / dL E \ = 

Sx l dx l dr \ dx' 1 J ’ 

together with appropriate initial conditions, can only determine trajectories in (Q n 
described by equations of the form 

Fxiq 1 , ■ • •, q n ) = 0, where K = 1,..., n — 1, (25) 

or, usually under obvious additional assumptions, of the form q K = q K (q n ). 

To determine the complete motion q l = q l [t) defined by the original Lagrangian 
L, one has to add to the n — 1 equations taken out from m the equation 

G(q i (t),tf(t))=E. (26) 

Thus, to determine the complete motion, one needs the triple (L e ,G,E). The pair 
(q l {t), t) geometrically represents a world line in the space of states Q ra x 1R in which 
the unit taken along the real axis IR is equal to the unit of the Newtonian time t. 
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Remark. Equation <P1) could as well be replaced by the equivalent equation 


<!>e (' q l {t),q>{t )) = 1 . 


5. The inverse Jacobi problem 

Let Lh (x 1 (t), x ' 3 (t)\ be a function homogeneous of degree one in the variables x ,% . 
A variational principle with Lh taken as the Lagrangian is only determining (non- 
parametrized) curves in a Q". The following questions can be asked here 


i. What data should be added to the knowledge of Lh, in order to be able to 
lift the spatial paths in Q™ to motions = q l (t ) determined by a Lagrangian 
L (q l (t ), q 3 (i)) such that the given homogeneous Lagrangian Lh is its Jacobi 
Lagrangian Le corresponding to E taken as the energy constant? 

ii. What is the algorithm that enables us to determine L in terms of an arbitrarily 
given Lh and what are the necessary additional data that make the solution to 
the problem unique? 


Problem of such a kind was formulated and solved in [ 5 ] under the name of inverse 
Jacobi problem. Now I would like to present its solution. 


All that said here so far suggests that a good candidate for the additional data 
would be an arbitrarily assigned function G(q l , q 3 ) being the hoped-for energy function 
of the yet unknown Lagrangian L. 


After introducing the velocity variable v l = q l (t ), relation (3.2) turns into a partial 
differential equation 

*■ ldL 1 'v n ^-L = G (27) 


dv l 


dv r 


for an unknown function L(v 1 ). In Eq. (12711 . G = G(v % ) is treated as a given function, 
and the dependence of L and G on q l is here suppressed. 

Applying the standard methods of integration of partial linear differential equa¬ 
tions, a general integral of m can be found to have the form 


L(q\v J ) = Vl g rs v r v s 1 l[ q\ V —— ,J\g P qV p v < i\\ +A(q\v 3 ), (28) 

\ vl 9kiv k v L | v ) 

where gij stands for the metric tensor in the manifold (Q ra (in case such a tensor is 
absent, one may write down = Sij), and where 

I(c\p):= J^Adp. (29) 


The function A (q z ,v 3 ) in (1281) is an arbitrary integration function homogeneous 
of degree one in the variables v 3 . The equation |28|) represents a general formula 
that determines a class of Lagrangians L describing a conservative dynamical system 
in terms of an a priori assigned energy function G of the system and an arbitrary 
homogeneous Lagrangian A. 
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To solve the problem, we have to remove the arbitrariness of A by making use 
of the requirement that the given homogeneous Lagrangian Lh(x l , x' 3 ) be the Jacobi 
Lagrangian corresponding to the Lagrangian L determined by Eq. P|). 

In order to be able to use the definition m of Le, we have to find first the 
function <j>E by solving the equation 

G (^x\ = E. (30) 


By using the requirement just mentioned, it is a quite simple technical matter to 
find the function A as a functional of L^, G, and 4>e- 

Substituting this functional into ca, one obtains the Lagrangian L which solves 
the problem posed: 


L{d\v j ) = yj 


gijV l vi 

-il q 


I q\ 


\/| g pq vPvi 

\A g rs v r i 


vl g rs v r v s 


V\ g P qV p v<i I ’ fan{q m ,v n ) 


+ L h (q i ,v j )-Ed> E (q i ,v j ). (31) 


6. Geodesics in a Lorentzian manifold 

Let g a /3 = <?a/?(£ 7 ), ot,(3 = 0, l,...,n, be a Lorentzian metric in a local coordinate 
system {£“} in a manifold M n+1 . The choice of its signature is + — . The 

geodesic lines £“ = £ a (f) in M" + , parametrized by an affine parameter f, are defined 
by the action 

T2 

W = J g a 0 U a u 13 dt, (32) 

Tl 

where u a = . The action (I32II determines geodesics as loci of points in an n + 2- 

dimensional space R x M” + , where R is the parameter axis. The space R x M Jl+1 is 
here, unlike in Newtonian mechanics, defined only locally over a geodesic line being 
just under consideration. In the case of the action (|32| . let us denote its “energy” 
function by G. By making use of Eq. m, the function G can be found in the form 

G(f \ u 13 ) = -\g a p u a vP. (33) 

If one assigns now to the “energy” constant C the value 

G = — \ em 2 c 2 , (34) 

where e = ±1, and m and c are some constants, then by solving the Euler-Lagrange 
equations, with m ^ 0, for e = 1 one obtains timelike, and for e = —1 spacelike 
geodesics. The assumption m = 0 is used here in case one wants to obtain a null 
geodesic. 

Since the Lagrangian in the action (|32|l does not depend explicitly on t , so it 
is possible here to perform the Jacobi reduction. To this end, one must first solve 
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Eq. (I2U1) in which G is replaced by G from Eq. lt221) . and E by C defined in Eq. (EH) . 
Thus the solution m takes now the form 

9\t) = (j) E {x a , x" 3 ) = ^ \J e g a p x' a x' 0 , (35) 

where x a {r) = £“ (9(t)) and x' a = Note that the Jacobi reduction is not possible 
in the case of null geodesics. 

Now with the aid of Eq. (1221) . the Jacobi Lagrangian Lq corresponding to the 
Lagrangian L(^ a ,x ,/3 ) of the action (13211 can be easily found as 

Lc{x a ,x ,/3 ) = —emc\jeg a ^x' a x'^. (36) 

The Lagrangian Lc{x a , x' 13 ) is homogeneous of degree one in x' a , with all the conse¬ 
quences of this fact which were indicated above. Thus, in case one would not like to 
introduce any additional constraint condition, geodesics can be described analytically 
only by equations of e.g. the form x l = x l (x°), i = 1 This means that the 

geodesics are loci of points in the manifold M" + , i.e., in this manifold, they are world 
lines in the terminology used in the theory of relativity. 


7. A Newtonian dynamical system 

Let us consider in (Q" a system defined by the Lagrangian 

C = \ e i3 uV + § A k v k - V ., (37) 

where the notation is a standard one; i,j = 1 It is assumed that the kinetic 

energy tensor eij, as well as the potentials Ak and V are functions of only the coordi¬ 
nates q l in Q", and they do not depend explicitly on the time t. The system satisfies 
then the energy conservation principle 

Q = \ uV + V = E . (38) 

If one wishes to apply to Eqs. m and the Jacobi reduction procedure described 
in Sec. 4, one has to solve with respect to 9' first the algebraic equation 


,/* 


x x 

W 


+ V = £ 


(39) 


corresponding in the present case to Eq. 113, and next to substitute into the equation 
which corresponds now to Eq. J23 the solution of Eq. (Ell) , which is 


9' = 4 >e{x k ,x' ) = 


2 {£ - V) ‘ 

The outcome of all the operations just described is the Jacobi Lagrangian 


(40) 


£s = sj2{£-V ) e^x' 1 ^ + ~ c A i x' i (41) 

of the system defined by the Lagrangian E3); for the notation cf. Eqs. mm- The 
Lagrangian m determines spatial paths in Q n , whereas the Lagrangian m is defin¬ 
ing motions in Q" which could be looked upon as world lines infix Q", where R is 
the Newtonian time axis. 
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8. Geodesics in a stationary space-time 

Let 


( 9 °) 

S = —smc 

( 9 °) 


1,9 J2 I - 

in 


dq k dq k dq l \ 0 

500 + 2 50 k d^ + 9 kl d^d^> dq 


(42) 


be an action for geodesics, q k = q k (q °), k = 1, n, in a space M n+1 with coordinates 
q a (a = 0,1, n). It is assumed that all g a p do not depend explicitly on q°. The 
minus sign is standing here to assure a principle of the least action, as well as the 
positive definiteness of energy. 

After replacing q° by t = —, the Lagrangian corresponding to m can be ex¬ 
pressed as 


L{q k ,v l ) = -me 2 e ye ^fifoo + 2 g ok — + g u ■ 

where v k = q k (t). Stationarity of M” +1 implies the conservation of energy 

goo + gok v k /c 


G = 


V £ (goo + 2g 0 kV k /c + g k i ( v k v l )/c 2 ) 


= E. 


(43) 


(44) 


Let us apply now the Jacobi procedure presented in Sec. 4 to the Lagrangian m 
taken together with its energy function m- It is a fairly straightforward matter 
to show that in case the function G is given by the expression (in . the algebraic 
equation m on O' has a unique solution of the form m in which for the function 
( f>E one must take 


where 


/ gok /k 

4>E = - X 

eg oo 


7 ij = 


E 


7 ^ x' l x'° 


CVS V E 2 / c 2 - m 2 c 2 e g 0 Q 

goigoj 


9ij 


goo 


(45) 


(46) 


is the so-called space metric tensor, cf. fX), and where the notation introduced in 
Eqs. and d applies. 

With the aid of Eq. Il22[l . the corresponding Jacobi Lagrangian Le can now be 
easily found as 


L e ( x\x' 3 ) = 


' E 2 
, c 2 goo 


m* c* £ l 7 ^ x' x 


Egok /fc 
egoo 


(47) 


The Lagrangian m determines spatial trajectories in M" being a section of M’ l+1 
with the hypersurface x° = const. 

Let us note that due to the equivalence principle, the mass parameter m that 
enters the Lagrangian BD , unlike the parameter e, does not appear in the Euler- 
Lagrange equations of motion which follow from this Lagrangian. These equations of 
motion admit however a whole class of solutions for which 


goo + 2 gok -h g k i 


v k v l 


= 0, 


(48) 
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i.e. for which L = 0. Thus these motions, represented by null geodesics in M n+1 , are 
not determined by an action principle based on the action ®. Solving the energy 
conservation law iTHli with respect to the square root of the expression standing on 
the l.h. side of Eq. ®, one can see that the square root tends to zero for to —> 0 for 
the values of e and E kept fixed and different from zero. Therefore, vanishing of the 
expression in I® can be considered to be equivalent to limm = 0. 1 

From Eq. it follows that Le is a meaningful Jacobi Lagrangian also for m = 0. 
Therefore, despite the fact that the action <E3> does not work for null geodesics in 
M" , the corresponding Jacobi action based on the Lagrangian Le, given by Eq. cm 
for to = 0, defines spatial paths in M" of such geodesics in JM" + . In that case Le 
is the Lagrangian of Fermat’s principle for stationary space-times. This principle, 
thought in a different theoretical framework, was already discussed e.g. in |5]. 

For to ^ 0, the action principle based on the Lagrangian Le given by m can 
be considered as being a generalization of Fermat’s principle for non-null geodesics 
in a stationary space time M" +1 . A Lagrangian of this kind, but only for static 
space times, was discussed in 0, which unfortunately is a paper with many logical 
and technical errors. In neither, however, of the papers just mentioned, the true 
dynamical origin of the principles discussed there was revealed. 

Let us finally observe that one can identify the manifolds M" and Q" of Sec. 7. 
This is implied by the fact that after making the identifications 


— Tij i 


eAi 

V 


90k E _ 

goo 

£ + b 


2 2 
TO C £ — 


E 2 


2 ffoo C 2 


(49) 

(50) 

(51) 


and fixing the values of to, e, £ , E, one can uniquely express the quantities g a p 
through eij, Ak, V, or vice versa. 

The identification of the spaces Q" and M" and the relations m-m demon¬ 
strate not only the equivalence of the two dynamics defined, correspondingly, by Le 
and C £ , but they also reveal the existence of maps leading from e.g. the dynamics 
determined by C to that by L or the other way round, in accordance with the diagrams 


A (A £) 


c £ , 3 c £ 

(G,E) 


(G , E) 


A 


(52) 


and 


L, (L, E) L e , 3 L e + (G, £) C, (53) 

(S, S) 

where the notation refers to objects that were already discussed in this article. 

To demonstrate the way of how diagrams of this kind should be read, let us explain 
it by taking the diagram 11521 ) as an example. The starting point here is the non¬ 
degenerate Lagrangian C of the form which describes the motion of a system as 

^■Of course, a similar equivalence could have been obtained by putting down m = 1 and letting 
e tend to zero. The way accepted in the article seems to be a more physical one. It demonstrates 
that a null geodesic is a limiting case of either timelike or spacelike geodesics which are selected by 
choosing either one of the two values of the discrete parameter e = d=l, while the mass parameter 
accepts its values from a continuous interval. 
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a locus of points in the n + 1-dimensional space R x (Q’\ The knowledge of C uniquely 
determines, by means of Eq. m, the energy function Q given by A choice that 
must be made before the Jacobi reduction procedure is started is that of selecting 
a value of the energy constant £ in Eq. m- Thus, to start the Jacobi reduction, 
one has to select a pair (£,£) in 1R x Q n . The outcome of the Jacobi procedure is a 
homogeneous Lagrangian Cs which is made equivalent to the Lagrangian Le by means 
of Eqs. m ~ m . Now, there exists a pair (G,E), consisting of a function G(x l ,x' J ) 
and a value of a constant E, such that when the piece of information encoded in the 
pair is logically added to that encoded in the Lagrangian Cg , one obtains the starting 
point of an inverse Jacobi procedure that leads us to the target Lagrangian L(q k ,v l ) 
given by Eq. m- Of course, for every Cs there is only one pair (G, E) that allows 
us to obtain the Lagrangian B31 1 . 

9. Geodesics in a stationary space-time in an affine parametrization 

The action C0 can be easily transformed to a homogeneous form. This may be 
achieved by introducing an additional dynamical variable q° (r) as a function of a new 
parameter r. Its values are here denoted by x°, i.e. x° = <jj°(t); and the remaining dy¬ 
namical variables are then transformed into x k = x k (r) := q k (< 7 °(t)) . After changing 
the integration variable q° —> t, dq° = x /0 dr, the integral (TTH takes the form 


r 2 

Sh = —emc J\Je g a /3 x' a x'^ dr , 


(54) 


where all g a /3 in the integrand do not depend explicitly on x°. The two actions, given 
respectively by m and IS}, determine the same loci of points in the space M n+1 
provided the function x°(r) in m is not a fixed one, but it is treated like any other 
dynamical variable during the variational procedure. This property of the action 
is due to the fact that its Lagrangian is a homogeneous function of degree one 
in the velocities x ,a . Thus the two dynamics, defined by the actions <E2> and |S} 
respectively, are mutually equivalent. 2 

There are various methods of transforming the action IS} into another one which 
would give us solutions in a form of parametrized curves in a suitably defined con¬ 
figuration space Q or, differently speaking, solutions which are world lines in locally 
defined spaces Ex Q, where ]R is the parameter axis. All such methods amount to 
adding new information to that encoded in the action IS}- For instance, one can 
substitute for the function x°(t ) in IS} any given, monotonous function x° = x°(r). 
This turns the action IS} into a non-homogeneous one which determines parametrized 
curves x k = £ k {r), where £ fe (r) := x k (x°(r)), in the space Q = M n . 

In this section, it is the inverse Jacobi procedure that is to be used. The Lagrangian 
of the homogeneous action IS} is 

Lh(x k ,x ,f3 ) = —emc^je g a px' a x'P, (55) 

2 In some old classical texts on differential geometry the action EJ is referred to as describing 
geodesics in an arbitrary parametrization. This phrase is, however, slightly confusing, for it is used 
to mean in a not yet specified parametrization. 
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i.e. it is formally equal to the Jacobi Lagrangian given by Eq. 13(ill . but now it does 
not depend explicitly on x°. The formal equality may indicate that the procedure 
we are going to use is, in a sense, inverse to the reduction procedure discussed in 
Sec. 6. Thus, seemingly, to trace the inverse procedure, it would be sufficient to read 
the equations of Sec. 6 in a reverse order from Eq. m to 03 . Although one could 
in this manner obtain a piece of helpful information, yet the inverse Jacobi method 
is more than this. It is in a way a procedure of lifting dynamics of a certain type 
from a configuration space Q to dynamics of a different type in the space K x Q, 
where 1R is the axis of a meanwhile unknown parameter t = 6(t ), which is implicitly 
introduced by a choice of an “energy” function. In principle, the choice of such a 
function is fairly arbitrary. In practice, however, this choice may be a guess based on 
the Jacobi reduction method applied to certain hoped-for target Lagrangians of the 
inverse procedure. 

In the present case, the starting homogeneous Lagrangian is Lh given by Eq. (1551) . 
and the space Q = M" + . In accordance with Sec. 6, the “energy” function is chosen 
to be 

G(£ fc V) := -±g a0 u a uP, (56) 

where g a g = g a p(£, k ), £“(i) = x a ((0 -1 (f)), u a = (§-, and t = 0(r), for 6'(r) ^ 0, is a 
new parameter. Also the choice of the “energy” constant C is, in principle, arbitrary. 
In order, however, to obtain a desired target Lagrangian, it is chosen, in accordance 
with Eq. G3, as C = — i em 2 c 2 . The next step consists in solving Eq. 03 adapted 
to the present notation. Its solution is presented in Eq. (1551) . After changing in the 
expression for (j)c{x k , x ,/3 ), given by Eq. (I551) . the names of the variables from ( x k , x ,/3 ) 
to (£ fc , u f3 ), we substitute this expression and that for the function G(£ fe , w/ 3 ) given by 
Eq. (1551) into Eq. 03) , to obtain the target Lagrangian in the form 

L{ey) = -\g a gu a u^. (57) 

i.e. a Lagrangian of the same form as that in the action 03 , but now the Lagrangian 
03 does not depend explicitly on £°. The parameter t introduced by the choice of 
the “energy” function 03 is an affine one. The Lagrangian 03 determines world 
lines in 1R x M" + , where 1R stands for the t axis. 

The Lagrangian 03 depends neither on t nor on £°. Its independence of t gave 
rise to the possibility of the Jacobi reduction procedure which was performed in Sec. 6, 
and here it would restore the starting Lagrangian Lh- Although the Lagrangian 03 
is independent of £°, it does depend on it 0 = so £° is a typical cyclic variable, and 

the existence of such a variable enables us to apply the Routh reduction procedure to 

the Lagrangian L as well. 

For the sake of this procedure, we replace in the Lagrangian L the names of the 
variables u a with v a and write down Eq. 03 in a way that explicitly exposes the 
dependence of L on the variable v°: 

L({ k ,vP) = -l (goo(v 0 ) 2 + 2g 0k v°v k +g kl v k v l ), (58) 

In order to eliminate from the Lagrangian (1581 the variables (x°,u°), we have to 
compute the quantity Vq defined in Eq. 0 for the case considered now. We have 


dl 

dv° 


V 0 := -goov 0 - g ok v k = p 0l 


(59) 
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and the solution for v° of the last equation above is 

v° = <f>(po,€(t),v J (t)) := - ^L. (60) 

goo goo 

The Lagrangian CD equals the Routh function 1Z (£*, cj>(po, £ fe , v 1 ), v\t) and takes now 
the form 

£po(£ fe V) = \ lkiv k v l +p 0 — v k + Po . (61) 

goo 2 mg 00 

The Lagrangian above is of the same type as that defined by Eq. GD• Therefore, we 
can identify the configuration space of dynamics defined by the Lagrangian El with 
the configuration space 1R x Q" of the dynamics discussed in Sec. 7. Comparing in 
the two Lagrangians, £ and £ Po respectively, the coefficients at the same powers of 
v k , we obtain 


e ij 

— 7 ij ) 

eAi 

gok 

= PoC — 

5oo 

V 

Po 2 

= —--h const 


2<?oo 


(62) 

(63) 

(64) 


Comparing next the two sets of relations, represented respectively by Eqs. iTTiI-ir>Tl) 
and by Eqs. El-El. we see that Eqs. Cil and m are identical, and Eq. m and 
(El can be made identical by assuming that po c = —E. Then Eq. (ii>i turns into 


E 2 

2 Soo c 2 


+ const, 


(65) 


which demonstrates that the two dynamics, defined respectively by C and £ Po , are 
fully equivalent. 

The content of this section is a generalization of a result by Eisenhart El who has 
shown that the trajectories of a general holonomic conservative system of n degrees of 
freedom in classical dynamics can be put into correspondence with geodesics of a suit¬ 
able Riemannian manifold S, where dim S = n+ 1. In 0, however, no use of methods 
of analytical dynamics was made, in particular of those concerning cyclic variables, 
but instead only tedious transformations of the underlying ODE were performed. 

Another reason which enables us to consider the result just obtained as a more 
general one than that of Eisenhart is that it permits one to prolong the sequence of 
mappings shown in the diagram El by a new sequence presented in the following 
diagram 


L tri ^ n L h , 3 L h + (G, C ) L, (L,p 0 ) C Po = £. (66) 

( G,C) 

The complete sequence, made by joining the sequences © and El one after the 
other, forms a closed loop. In an analogous way, with the help of the algorithms 
presented in this article, one can also prove the validity of the following sequence of 
mappings 

inverse 

£,3 £ + (Vo, po) L, (L, C ) L h = L. (67) 

CPo> po) 
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The composition 1 (15311 . (RTTll l of the corresponding sequences forms again a loop of 
mappings which passes exactly through the same dynamics as the previous loop, but 
the other way round. 

Thus the two loops, {El. Ell and {El. EZH taken together, define a com¬ 
muting diagram of mappings between all the dynamics discussed above. In terms of 
pairs consisting of Lagrangians and spaces of states 3 of corresponding dynamics, the 
diagram may be shown as 

(L, R xM" +1 ) x —♦ (L h ,M n+1 ) = (L, JM n+1 ) 

I I (68) 

(£, RxQ") «— (Cs , Q”) = (L e , M”). 


In this diagram the names of the procedures which labelled the corresponding mapping 
arrows, as well as other details concerning the definitions of mappings are suppressed, 
but they may be easily recovered by means of the diagrams El. El, HU- and cu¬ 
lt is rather remarkable that the seemingly arbitrary constant in Eq. El can be 
easily determined. This follows from the fact that the two dynamics, (£, R x M” + ) 
and (£, R x Q"), are invariant under translations of respective parameters t in the 
two configuration spaces. And this fact was not exploited yet. The invariance induces 
in the space R x M n+1 the conservation law: G(£ fe ,u^) = —^ em 2 c 2 , in accordance 
with Eqs. El and El- In order to project this conservation law on the space Rx Q", 
one must replace in it the variables «“ with i>“, eliminate from it the variable u°, and 
make use of Eq. pi . replacing p 0 c by — E. After all this is done, one obtains 


tu2 

- \ J = \ e m 2 c 2 . (69) 

Jg oo c 

On the other hand, the same invariance of the dynamics (£, R x Q”) induces the 
conservation law (138[) . Upon making use of the relation El , and eliminating from 
Eq. El the potential V by means of Eq. El, one transforms Eq. El into 

\ 7 ijV l v J - - - = E - const. (70) 

2300 c 2 

Eliminating now the kinetic term from Eqs. El and El. one finds that 


const = E + — em 2 c 2 , 


(71) 


which shows that the relations El and El agree with each other. 

The last result indicates that the mappings from the diagram El preserve various 
features of the three types of geodesics, labelled by the values of e and m. 
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